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We discuss the role of channel coupling in the surface properties of an inter-nuclear potential for 
heavy-ion reactions. To this end, we analyze the experimental quasi-elastic cross sectioirs for the 
-I- 13 q _|_ I05,i06pj systems using the coupled-channels approach by including 

the vibrational excitations in the target nuclei. While earlier studies have reported a negligible 
influence of vibrational excitation on the surface diffuseness parameter for spherical systems, we find 
a significant effect for the C-fPd systems. Our systematic study also reveals influence of transfer 


couplings on the surface diffuseness parameter. 

PACS numbers: 25.70.Be, 24.10.Eq, 25.70.Hi 

I. INTRODUCTION 

Even though fusion cross sections reflect the dynami¬ 
cal behavior of the nuclei involved in the fusion [1], ap¬ 
parently it is incapable to get insights into the interac¬ 
tion mechanism between nuclei simply by measuring fu¬ 
sion cross sections. That is, theoretical calculations, such 
as coupled-channels calculations (either semi-classical or 
quantal), are indispensable in order to understand the 
fusion dynamics [2-5]. These calculations, in compar¬ 
ison with the experimental data, address the influence 
of coupling between the relative motion and the nuclear 
intrinsic degrees of freedom, and thus the interplay be¬ 
tween nuclear structure and reaction dynamics [6, 7]. As 
a result, the coupled-channels formalism has become a 
powerful and standard tool in order to interpret experi¬ 
mental fusion cross sections [4]. 

One of the most important ingredients of such calcu¬ 
lations is an inter-nuclear potential. Various types of 
potential, such as a double folding potential [8, 9] and 
a phenomenological Woods-Saxon potential [10], have 
been employed. Among these, the Woods-Saxon form, 
described by the depth parameter, Vb, the radius param¬ 
eter, ro, and the surface diffuseness parameter, oq, has 
gained popularity due to its simplicity and ability to re¬ 
produce many experimental results. While representing 
the Woods-Saxon potential, its depth and radius parame¬ 
ters can be mutually adjusted to reproduce the Coulomb 
barrier of the system once the surface diffuseness pa¬ 
rameter is fixed. Obtaining the quantitative measure for 
the surface diffuseness parameter is therefore of crucial 
importance, as the curvature of the Coulomb barrier is 
largely determined by the surface diffuseness parameter, 
which by definition reflects the surface property of the 


* gkaur.phy@gmail.com 


internuclear potential. 

Three different methods have been used in literature 
in order to extract a diffuseness parameter from experi¬ 
mentally measured data. The first method is to analyze 
elastic scattering cross sections either with an optical po¬ 
tential model or with a coupled-channels approach. This 
methods has established a value of surface diffuseness to 
be around 0.6 fm for most of the systems [11]. The second 
method is to use high precision fusion cross section. It has 
turned out that this method leads to a much larger value 
of surface diffuseness parameter, ranging from 0.75 to 1.5 
fm, than the first method [12, 13]. Recently, Hagino et 
al. have proposed the third method, which uses a quasi¬ 
elastic (QE) excitation function at large backward angles 
[14]. They have suggested that at energies well below 
the Coulomb barrier, a QE excitation function is sensi¬ 
tive mainly to the surface property of nuclear potential 
and can be used to obtain quantitatively the surface dif¬ 
fuseness parameter. This method leads to the diffuseness 
parameter for spherical systems to be around 0.6 fm [15], 
while somewhat larger values of the diffuseness param¬ 
eter have been obtained for systems involving deformed 
nuclei. 

The different values of surface diffuseness parameter, 
obtained for spherical and deformed systems, lead to a 
realization that not only obtaining the quantitative mea¬ 
sure of diffuseness but also understanding how the cou¬ 
pling affects the surface diffuseness is of importance. For 
this purpose, several studies have been carried out in 
aiming at understanding the effects of different couplings 
on the surface diffuseness parameter. Casques et al. 
[17] have performed both the single-channel (sc) and the 
coupled-channels (cc) calculations for ^^S-|-^°®Pb, ^®^Au, 
^®®W, and ^™Eu systems and concluded that the rota¬ 
tion coupling shows a significant influence on the diffuse¬ 
ness parameter value. See also Ref. [18]. Washiyama 
et al. [15] have analyzed the QE measurements for 
32,34s_|- 208 p|^ systems with the single-channel calcula- 
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tions. They did not obtain any enhancement in the value 
of surface diffuseness parameter, and thus a significant 
effect of vibrational coupling on the diffuseness param¬ 
eter was not observed. Studies on many other systems, 
such as 16 ,i7,18q_|_ 92]y[Q have also shown similar re¬ 
sults. Hence, it has generally been considered that the 
vibrational excitation has a marginal effect on the surface 
diffuseness parameter. 

The aim of this paper is to investigate further the effect 
of vibrational excitations on the surface diffuseness pa¬ 
rameter. As in the rotational excitations, one may expect 
that the effect becomes more significant when the vibra¬ 
tional excitation energy is small. Furthermore, we also 
investigate the effect of transfer coupling, which has been 
established to play an important role in the fusion dy¬ 
namics [19]. To the best of our knowledge, no study has 
reported concerning how the transfer coupling affects the 
diffuseness parameter. For these purposes, we re-analyze 
the experimental data for I2,i3(^_|_i05,i06p(^ systems re¬ 
ported in Ref. [20] using the coupled-channels approach. 

The paper is organized as follows. In Sec. II, we first 
discuss the characteristic feature of the i2,i3Q_|_i05,i06pj^ 
systems with respect to the vibrational and the transfer 
couplings. We then detail the procedure to extract the 
surface diffuseness parameter through the analyses of the 
quasi-elastic data. In Sec. Ill, we present the results 
of our analyses and discuss the role of vibrational and 
transfer couplings. Finally, we summarize the paper in 
Sec. IV. 


II. METHOD 
A. Systems 

In order to study the effect of vibrational coupling on 
a surface diffuseness parameter, we re-analyze the exper¬ 
imental quasi-elastic scattering for the i2,i3Q_|_i06,i06pj^ 
systems, for which the data have been available in Ref. 
[20]. While Ref. [20] focused on extracting the bar¬ 
rier distributions, our aim in this paper is to extract the 
surface diffuseness parameter using the deep-subbarrier 
data. We choose these systems as the target nuclei, 
I05,i06p(j^ exhibit low-energy vibrational excitations {E* 
Ki 0.5 MeV for the first excited state in i04,i06p^^^ ppg 
projectile, along with the chosen target, is an ideal 
candidate to study the effect of transfer coupling due to 
the following characteristics: 

1. Positive Q-values for one-neutron transfer channels 
(see Table I), which are necessary for a system to 
exhibit significant influences of a neutron transfer 
coupling [19, 21]. 

2. The existence of the weakly bound valence nucleon, 
which ensures the transfer of neutron. 

3. Vibrational couplings in the target nuclei, which 
are in general weaker as compared to rotational 


TABLE I. The Q-value for the neutron transfer channels for 
the i2,i3Q_|_i05,i06pj^ systems, given in units of MeV. Here 
the negative and positive signs correspond to the pick-up and 
stripping reactions, respectively. 


System 

(-In) 

(-2n) 

(-tin) 

(-t2n) 

iUbp^ 

-9.161 

-15.744 

-2.148 

-3.953 

12c -t ^°®Pd 

-12.185 

-16.082 

-4.615 

-3.532 

13c -t ^°®Pd 

-t4.615 

-7.571 

-tl.082 

-7.681 

13c -t ^°®Pd 

-tl.59 

-7.908 

-1.385 

-7.261 


couplings so that the effect of transfer may not be 
significantly masked. 

In our systematic study, we carry out the coupled- 
channels calculations by including the vibrational cou¬ 
plings in the target nuclei in order to understand the 
influence of the vibrational coupling on the surface dif¬ 
fuseness parameter. Moreover, since the 
systems have a negative Q-value for the neutron pick¬ 
up reactions (see Table I), a comparison with the -I- 
I05,i06pj systems will elucidate the role of transfer cou¬ 
pling in the surface diffuseness parameter. 

B. Procedure 

To perform a systematic study, the single-channel and 
coupled-channel calculations have been performed using 
a scattering version of the CCFULL program [22]. For the 
coupled-channels calculations, we have included the dou¬ 
ble quadrupole phonon excitations in the target nuclei 
in the harmonic oscillator limit. The deformation pa¬ 
rameter and the excitation energy for ^*^®Pd are given 
by /3=0.229 and if*=0.512 MeV, respectively [20]. For 
the ^°®Pd nucleus, we have followed Ref. [20] and have 
taken the average in the adjacent nuclei, that is, ^°^Pd 
and ^°®Pd, which leads to /3=0.219 and A* =0.534 MeV. 

The nuclear potential used in the calculations has a 
real and an imaginary components, both of which are 
assumed to have a Woods-Saxon form. The imaginary 
part simulates a compound nucleus formation. We have 
chosen the strength to be large enough so that the flux 
does not reflect inside the barrier once the barrier is over¬ 
come. In the calculations, we have used an imaginary 
potential with the depth parameter of 30 MeV, the ra¬ 
dius parameter of 1.0 fm, and the diffuseness parameter 
of 0.3 fm. This choice of parameters confines the imag¬ 
inary potential well inside the Coulomb barrier with a 
negligible strength in the surface region. As long as the 
imaginary potential is confined inside the Coulomb bar¬ 
rier with a large strength, the results are insensitive to 
the parameters of the imaginary part. For the real part 
of the nuclear potential, the potential depth Vq is fixed 
to be 185 MeV. The value of radius parameter tq is then 
adjusted for a particular value of the diffuseness param- 
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eter such that the Coulomb barrier height Vb for each 
system becomes the same as that for the Bass potential 
[23]. This is possible because the effect of variation in Vb 
and tq on the Coulomb barrier height compensates with 
each other in the surface region. That is, for a given value 
of diffuseness parameter, the results do not significantly 
depend upon the actual choice of Vb, as long as the same 
barrier height Vb is maintained. 

To ensure that the barrier height for the single-channel 
and the coupled-channels calculations corresponds to the 
same value, we have slightly readjusted the potential pa¬ 
rameters for the coupled-channels calculations by using 
the fusion cross sections at energies above the barrier. 
The following steps have been taken for this purpose: 

1. For a chosen and fixed value of the diffuseness pa¬ 
rameter ao, and with the depth Vb = 185 MeV, the 
value of the radius parameter tq is determined such 
that the Coulomb barrier energy Vb reproduces the 
Bass barrier. 

2. The fusion cross sections, crfus\ are calculated using 
the single-channel calculation. 

3. The full coupled-channels calculations are then per¬ 

formed to obtain the fusion cross sections in the 
presence of the channel couplings, In gen¬ 

eral, even at energies above the barrier, these fusion 
cross sections are different from due to the po¬ 
tential renormalization [24]. By slightly adjusting 
the radius parameter, tq, we match the fusion cross 
sections Cf at energies well above the barrier with 

This results in a set of potential parame¬ 
ters (Vb, To, ag) that reproduces the fusion barrier 
height by taking into account the couplings to the 
intrinsic states with the coupled-channels calcula¬ 
tions. 

4. For every set of nuclear potential parameters, both 
for the single-channel and the coupled-channels 
analyses, the quasi-elastic scattering cross sections 
are computed. For the single-channel calculation, 
the quasi-elastic cross section corresponds simply 
to the elastic scattering cross section. On the other 
hand, for the coupled-channels calculations, the 
quasi-elastic cross sections correspond to a sum of 
elastic and inelastic cross sections. 

In order to find the best fitted value of the diffuse¬ 
ness parameter, the chi square method has been uti¬ 
lized. To this end, the data with dcrqei/daR > 1, where 
daqei/dfi and da^/d^l are quasi-elastic and the Ruther¬ 
ford cross sections, respectively, have been excluded from 
the fitting procedures, even though they are included in 
the figures for completeness. The uncertainty in the opti¬ 
mum value of ao has been calculated using the following 
procedure. For the Xmin value corresponding to the best 
fit value of the diffuseness parameter, the quantity (Xmin 
+ Xmin /’'V®'® calculated, where n denotes the number 


of degrees of freedom. The intersection of this quantity 
with the x^ envelope gives the two values ag and a[]" 
defining the error in ag. 

For the -h los.ioep^ 13^ + los.ioep^ systems, 
the QE excitation functions have been measured at 165° 
in the laboratory frame as reported in Ref. [20]. To 
be consistent, all the calculations have been carried out 
at same scattering angle. In order to ensure that the 
calculations are properly scaled according to the avail¬ 
able data, the calculated ratio of the quasi-elastic to the 
Rutherford cross sections are analyzed and plotted as 
functions of the effective energy defined as [25, 26], 


Eeff 


2Ec.ni. 

I -I- cosec 


( 1 ) 


where E^.m. and 0c.m. are energy and scattering angle 
in the center of mass frame, respectively. This corrects 
for the “angle dependent” centrifugal effects by making 

Crqe(Eeff, 180°) Ri CTqe (Ec.m. , ^c.m.) • 


III. RESULTS AND DISCUSSION 
A. Single-channel calculations 


Let us now numerically extract the surface diffuse¬ 
ness parameter from the experimental data for the quasi¬ 
elastic scattering. We have first performed the single¬ 
channel calculations without including the inelastic ex¬ 
citations of the target nuclei. The value of ag has been 
extracted using the procedure explained in the previous 
section. After the chi-square fitting, the best fitted value 
of ag for the -I- ^°®Pd and -I- ^°®Pd systems have 
been found to be 0.80±0.04 fm and 0.94±0.07 fm, re¬ 
spectively. The quasi-elastic cross sections obtained with 
several values of the surface diffuseness parameter are 
shown in Figs. 1 (a) and 1 (b) for the + ^°®Pd and 
-I- ^°®Pd systems, respectively. The chi-square fit for 
the -I- ^°®Pd system is shown in Fig. 2. Similarly, for 
the -h ^°®Pd (shown in Fig. 1 (c)) and + ^o^Pd 
systems (shown in Fig. 1 (d)), the best fitted values of ag 
after minimizing the have been found to be 0.64±0.05 
fm and 0.76±0.04 fm, respectively. 

The optimum values of the surface diffuseness param¬ 
eter are summarized in Table II. We notice that these 
values are significantly larger than the “standard value” 
of around 0.6 fm (obtained from elastic scattering cross 
sections), which are in a similar situation as in systems 
with a deformed target. 


B. Coupled-channels calculations 

In order to investigate whether the large values of 
surface diffuseness parameter obtained with the single- 
channel calculations are due to the neglect of channel 
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FIG. 1. (Color online) Comparisons of the single-channel calculations for the quasi-elastic excitation function obtained with 
several values of the surface diffuseness parameter, ag, in the nuclear potential. The panels (a), (b), (c), and (d) are for the 
- 1 - 12 q io®Pd, + ^°®Pd and -I- ^°®Pd systems, respectively. The experimental data are taken from Ref. [20]. 


TABLE II. The optimum value of the surface diffuseness 
parameter, ag, obtained with the single-channel and the 
coupled-channels calculations. Those values are given in units 
of fm. 


System 

single-channel 

coupled-channels 

r^C + iiJ^pd 

0.80T0.04 

0.69T0.04 

12c -b ^°®Pd 

0.94T0.07 

0.78T0.05 

13c -b “®Pd 

0.64T0.05 

0.60T0.05 

13c -b ^°®Pd 

0.76T0.04 

0.68T0.03 


coupling effects, we have next performed the conpled- 
channels calculations including the vibrational excita¬ 
tions in the target nuclei, los.ioepj^ value of oq 

has been varied and the best fitted value has been ob¬ 
tained after the minimization. The comparison of the 
coupled-channels calculations with several values of og 
with the experimental data is shown in Fig. 3. The best 
fitted value of og for the -I- ^°®Pd (shown in Fig. 
3 (a)) and the -I- ^°®Pd systems (shown in Fig. 3 
(b)) have been found to be 0.69±0.04 fm and 0.78±0.05 



FIG. 2. The chi square (x^) for the -I- ^'’®Pd system as 
a function of the surface diffuseness parameter, ag. The open 
and filled symbols represent the results of the single-channel 
and coupled-channel calculations, respectively. 


fm, respectively. Similarly, the best fitted value of ag for 
the -b i°’^Pd (shown in Fig. 3 (c)) and the -b 
i06pd systems (shown in Fig. 3 (d)) are 0.60±0.05 fm 
and 0.68±0.03 fm, respectively. Those values are sum¬ 
marized in Table II. 
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FIG. 3. (Color online) Same as Fig. 1, but with the coupled-channels calculations. 


It is apparent from Table II that the diffuseness pa¬ 
rameter decreases in the coupled-channels calculations 
as compared to the single-channel calculations with inert 
target nuclei. This observation is similar to that observed 
in the case of rotational coupling [17]. Notice that ear¬ 
lier studies with systems such as -I- ^°®Pb did not 
show any influence of a vibrational excitation on the dif¬ 
fuseness [17]. In contrast, our calculations for the -I- 
I05,i06p(j i3p; _(_ I05,i06pj systems show a significant 

influence even though the target nuclei are spherical. 

In order to understand this difference in the role of 
channel coupling, we proceeded to find the parameters 
responsible for the influence of a vibrational excitation on 
the surface diffuseness parameter. We have first checked 
whether the deformation of the nucleus is responsible for 
this behavior. The table III shows the deformation pa¬ 
rameter P for all the spherical target nuclei, ^°®Pb, ®^Mo, 
and i05,i06pj^ pQj. comparison, the table also lists the 
deformation parameter for the deformed target nuclei, 
and ^™Er considered in Ref. [17], which are es¬ 
timated from the measured B{E2) value [29] with the 
radius parameter of rp = 1.2 fm. As is shown in the 
table, ^°®Pb and ®^Mo have the deformation parameter 
which is significantly smaller than that for deformed nu¬ 
clei. For those nuclei, the vibrational effect on the surface 


TABLE III. The deformation parameter, f), the excitation 
energy, E*, and the spin-parity for the vibrational state in 
spherical target nuclei involved in the coupled-channel cal¬ 
culations. The values for the odd-mass nucleus, ^°®Pd, are 
estimated by averaging those quantities for the neighbouring 
^^'^Pd and ^°®Pd nuclei. The values are taken from the refer¬ 
ences listed along with. The table also lists those quantities 
for the deformed nuclei studied in Ref. [17], that is, and 
i"°Er. 


Nucleus 

P 

A* (MeV) 

state 

Ref. 


0.111 

2.615 

3" 

[28] 

®=Mo 

0.140 

2.849 

3" 

[27] 

lOSpd 

0.219 

0.534 

“2+” 

[20] 

106pd 

0.229 

0.512 

2+ 

[20] 


0.226 

0.112 

2+ 

[20] 

iTOEr 

0.336 

0.0786 

2+ 

[20] 


diffuseness parameter has been found to be marginal. In 
contrast, the Pd isotopes have a comparably large value 
of P to the deformed nuclei, leading to a large channel 
coupling effect on the surface diffuseness parameter. 

Naturally, one can expect that the energy of the excited 
state correlates with the deformation parameter. That is, 


























6 


the larger the deformation parameter is, the smaller the 
excitation energy will be. The table III shows the en¬ 
ergy of the first excited state of various nuclei considered 
in this paper as well as those in literature. It can be 
observed from the table that those nuclei which exhibit 
a large channel coupling effect on surface diffuseness pa¬ 
rameter have a small excitation energy of the first excited 
state, i.e., E* <1 MeV. In contrast, the first vibrational 
state of ^°®Pb is at 2.615 MeV, and thus the channel 
coupling effect is much smaller. 

Evidently, it is both j3 and E* which are responsible 
for the influence of channel coupling effect on the surface 
diffuseness parameter. The nature of coupling scheme, 
that is, the rotational versus vibrational, is unimportant 
with respect to the influence on the diffuseness. 


C. Role of Transfer Coupling 

In the coupled-channels calculation shown in the pre¬ 
vious subsection, the couplings to the quadrupole vibra¬ 
tional states are considered. If these were the only dom¬ 
inant channels, one would expect that the extracted sur¬ 
face diffuseness parameters were similar among the sys¬ 
tems. 

Fig. 4 shows the diffuseness parameter extracted with 
the single-channel (the open symbols) and the coupled- 
channels (the filled symbols) calcidations as a function 
of the mass product of the projectile and the target nu¬ 
clei, AtAp. It can be observed from the figure that the 
diffuseness parameter is reduced as the projectile isotope 
is changed from to (with the same target iso¬ 
tope) or the target isotope is changed from ^°®Pd to 
105pd (-^vith the same projectile isotope). This gives us a 
hint that there could be an effect of channel coupling in¬ 
volved other than the collective quadrupole excitations. 
As we have discussed in Sec. II, the transfer channel is 
a promising candidate for this, since a large probability 
of transfer is expected for due to the presence of the 
valence neutron. 

In order to investigate the role of transfer couplings, 
we have plotted in Fig. 5 the optimum value of surface 
diffuseness parameter as a function of the Q-value for 
neutron transfer. Since for the present systems, the two- 
neutron (2n) transfer is a second step process, the most 
important transfer channel is a one-neutron (In) transfer 
which may not be a general case. Hence, we have plotted 
the surface diffuseness as a function of the Q-value for 
the In transfer channels. For the In transfer channels, 
we have considered the -|-ln and the —In channels for the 
and projectile nuclei, respectively. These would 
be most preferable transfer channels from the point of 
view of the transfer Q-value, although the -I- In channel 
may be equally important for the ^^C-|-^°^Pd system. It 
can be observed from the figure that, as a general trend, 
the surface diffuseness decreases as the transfer Q-value 
increases. This might indicate that the difference in the 
surface diffuseness parameter between the projectile 


and the projectile could be attributed to the influ¬ 
ence of neutron transfer coupling. It would be an in¬ 
triguing future work to confirm this conjecture by carry¬ 
ing out coupled-channels calculations including both the 
collective excitations and the neutron transfer channels, 
although it is beyond the scope of the present paper. 


IV. SUMMARY 


The value of surface diffuseness parameter in the in- 
ternuclear potential for the -I- 13q 

_l_ I05.i06p(j systems has been extracted from the mea¬ 
sured quasi-elastic scattering cross sections at a back¬ 
ward angle. To this end, both the single-channel calcu¬ 
lations and the coupled-channels calculations including 
the quadrupole vibrational excitations in the target nu¬ 
clei have been carried out. Even though the systems 
studied involve spherical nuclei, the comparison of the 
coupled-channels calculations with the single-channel cal¬ 
culations revealed the reduction in the values of the dif¬ 
fuseness parameter. A similar reduction had been re¬ 
ported earlier for deformed systems due to the rotational 
coupling. Evidently, the conclusion in Refs. [17-19] that 
the vibrational excitation has a marginal effect on the 
surface diffuseness parameter is not a general one, but 
instead, the effect becomes significant even for the vi¬ 
brational coupling when the coupling is strong enough. 
We have argued that the reduction in the extracted dif¬ 
fuseness parameter is due to the strong coupling to the 
low-lying collective states, and the nature of coupling is 
not important. That is, the reduction can be observed 
both for the rotational and the vibrational couplings as 
long as the coupling strength is large enough. Further¬ 
more, a discussion has been made in order to under¬ 
stand the influence of the transfer coupling on the sur- 
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FIG. 4. (Color online) The extracted surface diffuseness pa¬ 
rameter for the 12,13Q_|_i05,i06pj^ systems as a function of the 
product of masses of the projectile and target nuclei, AtAp. 
The open and filled symbols represent the results of the single¬ 
channel and the coupled-channel calculations, respectively. 
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Q-value (MeV) 


FIG. 5. (Color online) Same as Fig. 4, but as a function of 
the Q-value for the neutron transfer channels. The lines are 
to guide the eye. 


face diffuseness parameter. We have observed that the 
surface diffuseness parameter gets smaller as the transfer 
Q-value increases. This implies that the surface diffuse¬ 
ness tends to be smaller when the transfer coupling is 
stronger. It would be an interesting future work to con¬ 
firm whether this trend holds in other systems as well. 
For that purpose, it would be interesting also to perform 
coupled-channels calculations including the transfer de¬ 
gree of freedom and carry out systematic studies in order 
to clarify the interplay between the collective excitations 
and the transfer couplings. 
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